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Abstract. An exact calculation of the Higgs boson decays H —> q;q :i mediated 
by flavor changing neutral currents is presented in the context of the Standard 
Model. Using up-to-date experimental data, branching ratios of the order of 10 -7 , 
10 8 , 10 8 , and Hr 17 are found for the bs + sb, bd + db, sd + ds, and cu + uc 
decay modes, respectively. 


Keywords: Higgs decays, flavor physics, Standard Model. 

1. Introduction 

In the light of the observation of a Standard Model (SM)-like Higgs boson [HE], some 
rare or suppressed processes of this particle become relevant. The Higgs boson has a 
relatively small decay width, and this is a good reason to revisit some of its rare decays, 
as the corresponding branching ratios may be significantly enhanced. Naturally 
suppressed (one-loop) decays of this particle are the flavor changing neutral currents 
(FCNC) decays H —> qiqj, with H the Higgs boson and = bs , bd, sd, cu. In these 

processes the corresponding amplitudes crucially depend on the quotient m\/m^q, 
with nik the mass of the quark circulating around the loop. For m|/m^ -C 1, the 
corresponding decays are strongly suppressed; this is known as the Glashow-Illiopolus- 
Maiani (GIM) mechanism [5J. Examples of such decays are the following FCNC 
transitions of the top quark: t —> cV, with V = j, g, Z (see SI), and t —1 cH (see [5]). 
However, no GIM suppression is present for rri^/rn^ > 1; this is the case for the well- 
known FCNC b —> sq decay, whose branching ratio is ten orders of magnitude higher of 
the aforementioned top quark decays. This suggests that the FCNC Higgs decays into 
quarks of down type, especially the bs mode, could have interesting branching ratios. 
Some features of these decays have previously been studied in the context of the SM. 
The qiqjH vertex was studied in [6] using some approximations. Several authors [7] 
explored the possibility of a very light Higgs boson via the b -A sH decay. Also, some 
technical aspects of the one-loop sdH vertex were studied in [8]. This problem has 
also been studied in some SM extensions: the SM with a fourth generation |S], the 
Two Higgs Doublet Model m and the Minimal Supersymmetric Standard Model nu. 
The purpose of this paper is to present exact formulae for the branching ratios of the 
FCNC H -a qiqq decays within the context of the SM and numerically analyze them 
with the use of up-to-date experimental data. 
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The paper has been organized as follows. In section [2j exact formulae for the 
H —> qiQj decays are derived in the context of the SM. Section [3] is devoted to discuss 
our results. Finally, in section [4] the conclusions are presented. 


2. The H —> qiQj decays 


In the unitary gauge, the H —> qjqi decay arises through the one-loop diagrams shown 
in figure [l] The corresponding invariant amplitude is given by 

3 

M(H qjqi) = y «(p»> Si) ( p L Pl + Fr Pr) v(Pj,Sj ), (1) 

where the projection operators are Pl,r = (1 T 75 )/ 2 , and the right loop amplitude 
Fpt is related to the left loop amplitude Fl via: 

Fr = F L (rm < — > mj). (2) 


Here 


Fl = \— E Wi 

V x w “ j 
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( 3 ) 


where V is the Cabibbo-Kobayashi-Maskawa (CKM) matrix, whereas Cq(1) and Bq(1) 
are Passarino-Veltman scalar functions given by 

Co(l) = C 0 {m 2 ,m 2 j,m 2 Hl , 


Co(2) = C 0 (m 2 ,m 2 , m 2 H ,m\,m 2 w ,rn%) , 

B 0 (l) = B 0 (m 2 ,rn 2 k) rnw ), 1 = 1,2; mi = mi, m 2 = m 3 , 
B a (r) = B Q (m 2 H ,m 2 ,m 2 ) , r = 3,4; m 3 = m w , m 4 = m k , 
B 0 { 5) = B 0 (f),m 2 k ,m 2 w ), 

B 0 (s) = Bq( 0, m 2 , m 2 ) , s = 6, 7; m 6 = m w ,m 7 = m k ■ 


( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 


The masses m, and mj are the external quarks’ masses, whereas m k denotes the mass 
of a quark circulating in the loop. 



Figure 1. Diagrams contributing to the H —> qiqj decay in the unitary gauge. 
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The correct implementation of the GIM mechanism requires that, before 
evaluating the loop amplitudes, those terms that do not depend on the internal quark 
mass must be removed, that is, terms of the form 


T. v ikVkj m Hl m w ) ; 


( 10 ) 


with ) being arbitrary functions which do not depend on the 

internal mass mk and which vanish for i j due to the unitarity of the V matrix. 
Taking into account this fact, we can write the various form factors appearing in the 
loop amplitudes as follows: 

x k (x k + x w - Xj) 
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— Xk \xi(xi — l ) 2 + XiXj (3 Xi — Xj ) + Xj (l — 3x 2 — 6 Xj) ] 

— Xj [xi (l — x 2 ) + XiXj (2 Xi — Xj + 2) + 2 Xj (xj — 1)] |, (15) 
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where we have introduced the dimensionless variables x,- 


(16) 

(17) 

(18) 

(19) 

( 20 ) 


= ™t/r 


= m 2 /r 


Hi '"'Hi 

2 In 


Xk = rnl/mjj, and Xw = 7n^ v /m 2 H . In addition, <5=1 — 2(xi + Xj) + (xi — Xj) 

/o, terms which do not depend on mk have been removed. We proceeded in the same 
fashion for the form factors g 3 and gg, as the associated Bq(3) and Bq(6) functions 
do not depend on nik■ It turns out that the loop amplitudes are free of ultraviolet 
divergences, as it can be verified by adding together the form factors associated with 
the Bq(1) functions. This leads to 
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\ ' Xj — Xj 

9i = 


1=1 


4 t 2 


( 21 ) 


which vanishes after using the unitarity of the CKM matrix. In this respect, the 
numerical evaluation of the Bq functions must be performed with some care, as these 
functions may still contain a mfc-independent part, which could lead to additional 
cancellations due to the GIM mechanism. This means that we cannot use software 
such as FF m or LoopTools Mi to directly perform such evaluation. Instead, we 
use their analytic solutions in order to remove any redundant contribution to the loop 
amplitudes. As it can be seen from Eqs. © to ©. they include three distinct B 0 
functions: (i) In Eq. ©, -Bo(l) and B 0 (2) have an analytic solution, namely 


B o (0 — A + 1 + 


Xw 


■ log(xw) - 


Xk 


■ log(xfc) + F(xi,Xk,xw) 


( 22 ) 


Xk — Xw Xk — Xw 

where Z = 1,2, Xjq 2 ) = x (i,j)i an d A = i — 7 ^ — log ^ 7772 ^ is a divergent factor, which 
is common to all the Bq functions and thus vanishes in the total amplitude, since this 
is free of ultraviolet divergences. In addition, the function F(xi,Xk,xw) is 


F(xi,Xk,xw) = 1 + ^ 


1 ( x k Xw 


Xk 


Xi 


x+x_ 


log 
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x+ + x 
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V Xk 
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Xi \x + - X_ ; 

where x± = \J (yfxf ± y/xw) 2 ~ Xi- Notice that F(xi,Xk,Xw) vanishes in the limit 
0. (ii) In Eq. ©, functions Bq(3) and Bq(A) also have an analytic solution which 


xi 

is 


B 0 (r) = A + 2 - log(x) - F(x), 


(24) 
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where r = 3,4, the variable x stands for xk or and 

F{x) = 2\/4x — 1 arctan ( ) , a: >1/4. 

V V4x — 

Finally, (iii) functions B 0 (5), B 0 (6) and B 0 ( 7) in Eqs.([E|) and © become 

B 0 ( 5) = A + 1 H-—— \og{x w ) -—-log(a: fe ), 

Xk — Xw Xk — Xw 

B 0 (6 ,7) = A - \og(x w , x k ) ■ 

Once plugging the expressions for the Bo functions into the Fl amplitude, 
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(26) 
(27) 


where 


+ h log 


( Xw \ 
\x k J 


i=i 

UF(x w ) + f 5 F{xk) 


+ f 6 F{xi,x k ,x w ) + f 7 F(x j: x k ,x w ) 


fo = - j [{Xk + X W )i 1 - Xi + Xj) + Xi - 1 + Xj{Xi - 


h = 

h = 
k = 


Xw <5 
Xk 


'J )\ 1 


[{x k - Xw){Xi - Xj - 1) +Xj{xj - Xi - 1)] 


2xw{xk — xw)h 
h = 2xw§ ( 2xw + ~ x o ~ !) ■ 


Xk 

2 xwh 
Xk 


[2 {xk + 2xw)0 — Xi + Xj) + 2Xj{xi — Xj) + Xi + xf] , 


Xw) 


AxiXwD 

x [xi (4xi(l — Xi) — XiXj + 2xj{3xj — 1)) — Xj( 1 — Xj) 2 ] 

+ 2xi \xj iixi + Xj) 2 — 2xj(2xj + 1) — l) + 4xj(l — Xj)] |, 


(28) 

(29) 

(30) 

(31) 

(32) 


(33) 


fr = 


Xk 


+ xw) [(1 — Xi — Xj) 2 + Axj{xi — 2Xj)] 

4 


4 xwD 
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(34) 


where D = {xi — Xj) 3 — 2{x 2 — x 2 ) + Xi — Xj. The Fl amplitude is now free 
of redundances. The right loop amplitude is obtained from Fr = Fl(xi <->■ Xj) 
[cf. Eq. ©]. 


3. Discussion 

It can be verified that our results satisfy the following consistency conditions: (*) The 
amplitudes for the diagrams (1) and (2) vanish for m; = 0 = rrij (diagrams (3) and 
(4) do not exist in this case); {ii) loop amplitudes Fl and Fr are free of ultraviolet 
divergences once the unitarity of the CKM matrix is used; ( iii ) in the limit rrij — > 0 
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(rnj —> 0) the right (left) loop amplitude F R (F L ) vanishes and the left (right) loop 
amplitude is different from zero. 

We now explore the behavior of the loop amplitudes in the heavy mass limit. 
Let m\v = vnk in the expression for F R given by Eq. ©, in this limit we have 
C'o(l) = C 0 (2) = C 0 (m?,m?,m%,rnl,rnl,rnl), B 0 ( 3) = B 0 ( 4) = B 0 {m 2 H , m|, ml), 
and B 0 (5) = B 0 (6) = B 0 (7) = B o (0,m1,ml). Therefore 


where 


In addition, 


F l 


my^—^rrik 


log(z fe ) 


Cl2f{Xk 


V 1 4 V X k 

k T 

+ biF(xi,Xk) + b 2 F(xj,Xk) + b 3i F(xk ) j , 


(35) 



b 2 = 


b 3 4 = 



(36) 


(37) 


(38) 

(39) 


f{Xk) 


F(x,x k ) 



X — X^Xj 2 


(40) 

(41) 


where f{xk) is the solution of the function for 

nii = 0 = mj, provided that this function depends weakly on these masses. From 
these expressions, it is easy to show that each term within FL\m w ^m k separately 
vanishes for Xk —> oo. 

Let us now provide the calculations for the branching ratio of the H —> qjqi + qiqj 
decays. In general, one has 


BR(H —>■ q^qi + qtqj) 



(1 - Xi- 


xj)(\F l \ 2 + \F r \ 2 ) 


- Ay/XiXjRe (F l F r ) , 


(42) 


where Nc = 3 is the color index and we have added a factor of 2 in order to consider 
both, qjqi and qiqj, possibilities. To evaluate all the branching ratios of the allowed 
decays, the following values for the Higgs boson mass and decay width were used: 
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run = 125 GeV and Fr = 4.403 x 10” 3 GeV; the values for the remaining parameters 
involved are those reported by the Particle Data Group |T3] . 

We evaluated the three point scalar functions Co(l) and Co(2) using FF [T2] 
and LoopTools [T3] . Also, we perform an independent numerical evaluation of 
these functions by starting from their integral representation. In terms of Feynman 
parameters, the C'o(l) function can be written in the following form: 


m^Co(l) = f dx 

Jo 


2 

7KL 


arctan 


(2 — a)x — 1 

71 


— arctan 


1 — ax 


(43) 


where a = 1 + — Xj and 


A = 4 xw — 1 + 2 (1 — Xi — Xj + 2(xk — xw)) % — 6 : 


(44) 

Some simplifications are obtained in the limit Xi = 0 = Xj. In this case, the above 
integral reads as follows: 


'o(l) — - / 
Jo 


dx 


7a 


arctan 


(1 — x 

, 7a . 


(45) 


where 


A = 4 xw — 1 + 2 (1 + 2(xk — xw)) x — x 2 


(46) 


However, as it occurs in the exact case, this integral cannot be expressed in terms of 
elementary functions. If in addition, one assumes that xw = %k, the result (TfOl) is 
obtained. On the other hand, the corresponding expression for rnjjCo(2) is obtained 
from m^Co(l) via the interchange Xk O xw■ The numerical evaluation of the integral 
given by Eq. (TTH) leads to results that are in excellent agreement with those obtained 
using the FF and LoopTools programs. 

The branching ratios have been evaluated using the following approximations. 
For those decays of type down, only the contribution of the quark top was considered 
for the channels hs and bd , whereas for the ds channel, only the contribution of the 
c quark was taken into account. It results that in the case of the ds channel, the 
t contribution is quite marginal with respect to the c contribution due to a strong 
suppression factor coming from the CKM matrix. While the function Fl {Fr) induced 
by the t quark is one order of magnitude (of the same order of magnitude) with respect 
to the c contributions, the CKM effects are, respectively, \V st V t d\ = 3.36 x 10 -4 and 
IVgcWdl = 0 . 222 . As far as the cu channel is concerned, only the contribution of the 
b quark was included. 

Our results are displayed in Tabled] The values of itirC 0 (1) and m 2 H C 0 (2) are 
shown in Tabled] From Tabled] we can see that the FCNC Higgs decays into the bs , 
bd, sd, and cu modes have corresponding branching ratios of 3 x 10 -7 , 1.14 x 10 -8 , 
1.2 x ICC 8 , and 5 x 10 -15 . It is worth comparing these branching ratios with those 
associated with the decays b —> sy, t —> eg, t —> cy, t cZ, and t —> cH, which are 
of the order of 10 -4 , 4 x 10 -11 , 5 x 10 -13 , 10” 13 , and 10 -14 , respectively. Although 
significant, compared with the FCNC top quark transitions, the U —>■ bs decay is quite 
suppressed to be detected in future experiments. 

Our results should be compared with recent constraints for flavor violating Higgs 
decays derived from low-energy data. By assuming a renormalizable general effective 
Lagrangian for the Yukawa sector, bounds on lepton and quark flavor violating 
decays of the Higgs boson were analyzed in references [15] dS]. In particular, 
experimental limits on B s physics were used to derive constraints on the Fl —» bs 
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Table 1. Branching ratios for the FCNC Higgs decays H — » q t qj. No 
approximations were made. 


H -A q^j 

\V ik V kj [* 

F l 

Fr 

BR 

H -A bs 

1.66 x 10” 1 2 3 4 

-0.2246 

0.1769 

2.92 x 10- y 

H -A bd 

7.34 x 10" 5 6 7 

-0.2247 

0.1469 

1.14 x 10" 8 

H -A sd 

4.9 x 10~ 2 

7.95 x 10" 4 

-1.06 x 10 -2 

1.2 x 10" s 

H -A cu 

2.88 x 1CT 8 

(-3.8 + 0.047*) x 10" 3 

-8.41 x 10 -3 

5.31 x 10~ 15 


Table 2. Values of the diverse scalar Cp functions. 


H —> qiqj 

m-H Co(l) 

mjj Co (2) 

H -+ bs 

-0.79916758 

-0.40298304 

H -+ bd 

-0.79916754 

-0.40298302 

F[ -+ sd 

-3.17017656 

-0.58 - 3.86* 

H -+ cu 

-3.13929925 

-0.61-3.85* 


decay. The authors of these papers found that BR(H —> bs) < 4 x 10 -4 [TB] and 
BR(H —> bs) < 2 x 10 -3 [16] . Although away from the SM prediction found here, 
these branching ratios would even be out of the reach of the LHC due to large QCD 
background. 

4. Conclusions 

In conclusion, we have presented exact formulae for the FCNC Higgs decays H -+ qiqj 
in the context of the SM. Recent experimental data were used to predict the branching 
ratios for all the kinematic allowed modes. Although the branching ratios for the decay 
H -A bs are significantly larger than the ones associated with the top quark transitions 
t —> Uij, t —> UiZ , and t —> UiH , our numerical analysis suggests that these decays are 
out of reach of future experiments, and thus they may be very sensitive to new physics 
effects. Recent analysis using experimental data from B s physics shows that branching 
ratios up to four orders of magnitude larger than the SM prediction could be allowed 
for BR(H — > bs ); however, they would still be undetectable at the LHC [13 HI]. We 
think that our results can be useful for people interested in investigating these decays 
in other contexts of new physics. 
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